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Dark energy perturbations are normally either neglected or else included in a purely numerical 
way, obscuring their dependence on underlying parameters like the equation of state or the sound 
speed. However, while many different explanations for the dark energy can have the same equation of 
state, they usually differ in their perturbations so that these provide a fingerprint for distinguishing 
between different models with the same equation of state. In this paper we derive simple yet accurate 
approximations that are able to characterize a specific class of models (encompassing most scalar- 
field models) which is often generically called "dark energy" . We then use the approximate solutions 
to look at the impact of the dark energy perturbations on the dark matter power spectrum and on 
the integrated Sachs- Wolfe effect in the cosmic microwave background radiation. 

PACS numbers: 98.80.-k; 95.36.-fx 
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I. INTRODUCTION 

More than ten years after the supernova observations 
P, led to the general acceptance that the expansion of 
the Universe is accelerating, we are still far from a consen- 
sus on what is responsible for the acceleration. Although 
there is a name for the phenomenon, dark energy, there 
is as of yet no convincing physical explanation. 

For this reason it is imperative to learn as much as 
possible from observations. The fundamental observable 
considered so far is the equation of state parameter w 
which connects the average energy density p(t) and the 
average pressure p(t) through p(t) = w(t)p(t). This pa- 
rameter characterizes the background expansion rate and 
the distances. While some models make specific predic- 
tions, like iv = — 1 for a cosmological constant, there are 
whole families of models that can lead to any desired evo- 
lution of w(t) (possibly with some weak constraints like 
w > —1). Typical examples include scalar-field mod- 
els like Quintessence or K-essence, or generalized gravity 
models like scalar-tensor and f(R) theories. These the- 
ories cannot be ruled out based on a measurement of 
w alone. However, in general their perturbations evolve 
differently @, H i, @, 0, S, Q . The dark energy pertur- 
bations are therefore like a fingerprint of these models 
and, if measured, will allow to discriminate much more 
precisely between competing theories, and hopefully will 
allow to shed some light on the physical nature of what- 
ever accelerates the expansion of the Universe. 

For many models, the behavior of perturbations at 
the linear level can be described in terms of those of 
a fluid with a certain sound speed. This is the case 
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for Quintessence (canonical scalar-field) models for which 
the sound speed is c 2 s = 1, and for many K-essence models 
where the sound speed is arbitrary. In this paper we con- 
centrate on this class of models, and assume in addition 
that both Cg and w are constant. This latter assumption 
is usually violated, but the results should nonetheless al- 
low insight into the behavior of the perturbations. For 
models where the quantities vary only slowly with time, 
we expect the results to still hold in an averaged sense, 
due to the indirect nature of most observations. The big 
advantage of making these assumptions is that it allows 
us to solve the perturbations analytically under the ad- 
ditional condition of matter domination, leading to sur- 
prisingly simple results. In addition, when expressed in 
terms of the change of the gravitational potential rela- 
tive to the case without dark energy perturbations then 
the simple formulae turn out to be a surprisingly good 
approximation until today. Our results should be seen 
in this context, as it is of course always possible to solve 
the perturbation equations numerically. However, ana- 
lytical results allow a much better insight and also an 
easy way to see how the behavior changes as a function 
of the parameters. 

Experiments tell us that the dark energy is at most 
very weakly coupled to the things that we can observe 
directly, like galaxies and the cosmic microwave back- 
ground (CMB). It is therefore not only necessary to de- 
termine the dark energy perturbations but also to con- 
nect them to actual observables. We use our simple an- 
alytical results and make a small step in that direction, 
trying to establish the impact of the dark energy pertur- 
bations onto the matter power spectrum and the inte- 
grated Sachs- Wolfe (ISW) effect in the CMB. Our treat- 
ment here is far from complete and only an initial at- 
tempt to deal with the observational impact of the dark 
energy perturbations. In a follow-up publication [ic| we 
will pay particular attention to the question whether it 
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is possible to observe the class of perturbations that we 
study here, using probes of lensing and of the matter 
power spectrum. 

In detail, the paper is organized as follows. We begin 
with a short discussion of the perturbation equations to 
set the scene and to define our variables, and remind the 
reader of the solution for the matter perturbations during 
matter domination. We then derive simplified solutions 
for the dark energy perturbations during matter domina- 
tion. Those expressions are a good fit in their respective 
domain of validity, but once the dark energy starts to 
dominate they deviate from the numerical solution. We 
then show that the function Q(k, t) used in |j is well 
described by our solutions, even after matter domination 
ends. This is the main result of the paper. We finally 
use this observation to consider the impact of the dark 
energy perturbations on the matter power spectrum and 
the growth rate of the matter perturbations, as a func- 
tion of w and c 2 . We also investigate which aspect of the 
dark energy perturbations affects the ISW effect most 
strongly, before concluding. The appendices give more 
details on how we compared our analytical formulae to 
numerical results and on the evolution of the decaying 
modes that were neglected. 



II. FIRST ORDER PERTURBATIONS IN 
MATTER AND DARK ENERGY 

Throughout this paper, we will use overdots to denote 
derivatives with respect to conformal time r, related to 
cosmic time by dt = adr. We will denote the physi- 
cal Hubble parameter with H and with 7i the conformal 
Hubble parameter. We consider only spatially flat uni- 
verses, and our metric convention is defined by the line 
element, 



ds 2 



(1 + 2ip) dr 2 + (1 - 2(f) dx t dx l 



(1) 



We are therefore working in the Newtonian or longi- 
tudinal gauge, which influences the resulting perturba- 
tions especially on scales larger than the Hubble horizon, 
k<aH. On much smaller scales the choice of gauge is 
less important, and observables are independent of the 
gauge choice. 

The perturbation equations for a fluid with equation 
of state parameter w = p/p are 



5' = 3(1 +w)(j)' 



Ha 2 a ( p 



(2) 



V = -(l-3w)- 
a 

-(1 + -)^- 



k 2 5p . k 2 , 

p + (l + w ) V (3) 



Ha 2 p 



Ha 2 



where S — 5p/p is the density contrast, V = ikjTy/p is 
the scalar velocity perturbation (see also [l3j]) and the 



prime means the derivative with respect to the scale fac- 
tor a. In this paper we will look only at fluids with van- 
ishing anisotropic stress, a — 0, so that there is a single 
gravitational potential, 



3aH 



(4) 



and ip = <f>. The sum on the right hand side runs over all 
fluids. For our purposes, we will assume the presence of 
a matter fluid with w = Sp = as well as a dark energy 
fluid, parametrised by a constant w and a sound speed 
c 2 which determines the pressure perturbation through 



Sp = c s pS 



3aH i 



k 2 



-pV. 



(5) 



where c 2 = p/p is called the adiabatic sound speed of 
the fluid. This is not always a good parametrization, for 
example for models that cross the phantom barrier w = 
— 1 [l3[ or when mimicking modified gravity models Q 
(in which case also generically a ^ 0). However, it covers 
a wide class of models, for example canonical scalar fields 
(Quintessence, in which case c 2 — 1) and other scalar- 
field models (like K-essence, which allows for c 2 ^ 1). As 
we only consider models with constant w, we have that 
= w. The perturbation equations @ 



w 



3H(l+w) 

and ((4| become in this case 



6' 



V 



V 
Hm 2 



9a 2 H 2 (c 



k 2 



- (c 2 s -w) 6 + 3(1 + w)<f>' 



V k 2 c 



2^,2 



Ha 2 



(1 



Ha 



(6) 
r0 (7) 



In general it is difficult to solve these equations. To 
simplify the problem, we will assume that the universe is 
matter dominated. This means that the Hubble expan- 
sion rate is given by 



H 2 



H Q, 



8ttG 



-Pn 



(8) 



and that only the matter perturbations contribute to the 
gravitational potential in Eq. ((4]). We can therefore first 
solve for the matter perturbations alone, without consid- 
ering the dark energy, and then use the resulting gravi- 
tational potential as an external source in the equations 
for the dark energy perturbations. Of course this ap- 
proximation, and even more the assumption of a matter 
dominated background evolution, will change the results, 
and we will need to study what happens as the assump- 
tions break down. 

The perturbation equations for the matter perturba- 
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tions then are 

Kn = 

V = 

v rn 

k 2 <t> = 



Vm_ 

Ha 2 

Vrn 
a 

2a 



-3<f>' 

k 2 
Ha 2 ' 

5 m 4 



3aH 



-V„ 



As is well known, a solution to this set of equations is 



S m = S n ( a + 3 
V m 

k 2 6 : 



J L rn 
k 2 

-6 H Vn^a 1/2 



= 5 n a 1 



H 2 a 2 
' k 2 





10,0 


(9) 
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(10) 
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(12) 





(13) 
(14) 



where the constant S sets the overall scale (since the 
equations are linear). This can be verified simply by in- 
serting them into the differential equations. The value 
of 5 is set by the initial conditions and is in general 
a function of k. Additionally the growth of the matter 
perturbations is suppressed during radiation domination, 
so that S (k) is smaller than expected for the scales that 
enter the horizon during radiation domination. However, 
the same happens for the dark energy perturbations so 
that the effects cancel out when comparing the pertur- 
bations in matter and dark energy. For more details see 
Appendix \X\ which discusses how we compare the ana- 
lytical results to CAMB [lH ]. 

We see that the gravitational potential is constant as 
a function of time on all scales. The matter perturba- 
tions grow linearly with a on small scales, and they are 
constant on super-horizon scales (but their behavior on 
those large scales depends on the gauge choice). 



III. SOLUTIONS FOR THE DARK ENERGY 
PERTURBATIONS DURING MATTER 
DOMINATION 



We will now use the constant k 2 cf> of the last section to 
look for solutions to the general perturbation equations. 
We will study them in different limits, and then compare 

See also [H| 1 



for 



the results with full numerical solutions 
results using different approximations. 

Generically, we expect at least three regimes with dif- 
ferent behavior of the perturbations: 

• Perturbations larger than the causal horizon 

• Perturbations smaller than the causal horizon, but 
larger than the sound horizon 

• Perturbations smaller than the sound horizon 

We start by looking at perturbations larger than the 
sound horizon, k <C aH/c s . In this case, we neglect all 
terms containing the sound speed in Eq. ([7]), effectively 




FIG. 1: The figure shows the behavior of the variable $de- 
The black dot-dashed line is the numerical solution with 
c 2 3 =0.1 and w = -0.8 for the mode k = 200H . The red 
solid line is the approximation on scales above the sound hori- 
zon, Eq. (|16|l and the blue dashed line is the approximation 
to the scales below the sound horizon, Eq. fTTJ). The two ver- 
tical lines give the scale factor at which the mode enters the 
Hubble horizon (left line) and the sound horizon (right line). 
The numerical solution shows how the perturbations decay at 
late times when matter domination ends, but radiation was 
omitted from the numerical calculation to allow for a longer 
dynamic range in a to illustrate the different regimes. 



setting c 2 — 0. The solution for the velocity perturbation 
is (neglecting a decaying solution ocl/a) 



V = -S (l+w)H 



(15) 



Up to the prefactor (1 + w) this is the same as for the 
matter velocity perturbations. We find that this expres- 
sion is valid on scales larger than the sound horizon even 
if the sound speed is non-zero. 

We can now insert this solution for the dark energy 
velocity perturbation into Eq. ©. Again setting c 2 — 
we find the solution 



5 = 5 (l + w) 



a 



1 - 3k; 



(16) 



where we neglected a term proportional to a 3w which is 
decaying as long as w is negative. Not surprisingly, also 
this solution becomes equal to the one for matter pertur- 
bations for w — > 0. Relative to the matter perturbations 
the dark energy perturbations are suppressed by the fac- 
tor (1 + w). This factor is necessarily always there, as 
the gravitational potential terms contain it. It can be 
thought of as modulating the strength of the coupling of 
the dark energy perturbations to the perturbations in the 
metric. For w — —1 the dark energy perturbations are 
completely decoupled (in the sense that they do not feel 
metric perturbations - but they can still produce them 
if the dark energy perturbations are not zero). 

Just like Eq. (fT2")) the solution is composed of two 
terms, where the second one dominates on scales larger 
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than the horizon (k < aH) and the first on smaller scales 
(but still larger than the sound horizon). The super- 
horizon part of the solution is the same as for matter 
apart from the overall (1 + w) factor, while the sub- 
horizon solution is additionally suppressed by a factor 
of 1/(1 — 3w) relative to the matter perturbations. 

Although these factors can suppress the dark energy 
perturbations significantly compared to the dark matter 
perturbations, especially if w is close to —1, the existence 
of a sound horizon is even more important. Inside the 
causal horizon, the dark matter perturbations grow lin- 
early with a (until the perturbations become non-linear). 
The dark energy perturbations on the other hand will 
eventually encounter their sound horizon if c 2 > 0. Once 
inside the sound horizon, they will stop growing. This 
means that the dark energy perturbation spectrum is cut 
off on small scales. 

To get a solution on small scales, k ^> aH/c s , we start 
again with the equation for the velocity perturbation. 
However, we expect the two terms with k 2 to cancel to 
a high degree to avoid large velocity perturbations, or in 
other words 



(l + w)(f> 



H 2 Q, m 



(17) 



We find that the dark energy perturbations stop growing 
and become constant inside the sound horizon. 

The velocity perturbations are now given simply by 
using Eq. Q and inserting Eq. (fT7|) : 



V 



-3Ha(c 2 s - w )S = - 9 -(l + w){c 2 s W )^^a-^ 2 . 

(18) 

is not important 



The extra term in brackets in Eq. 
for the scales of interest here. 

Finally, we would like to remind the reader that these 
results have been obtained under the assumption of a 
time- independent w and c 2 . On the other hand, a k 
dependence of c 2 is allowed. 

As the horizons grow over time, a fixed wave num- 
ber k will correspond to a scale that is larger than the 
causal horizon, k < aH, at early times, and eventually 
it will enter the causal horizon and later the sound hori- 
zon. This makes it possible to illustrate the behavior 
of the perturbations in the different regimes in a single 
figure: In Fig. [T] we plot the numerical solution for the 
dark energy density contrast for k = 200iJo as well as 
the expressions (TIB)) and (TIT)) . It is easy to see how the 
perturbations start to grow inside the causal horizon but 
how the growth stops when the sound horizon is encoun- 
tered and pressure support counteracts the gravitational 
collapse. 



IV. DARK ENERGY DOMINATION AND Q 

Matter domination was a crucial ingredient to compute 
the behavior of the dark energy perturbations, since the 



gravitational potential <f> is constant while matter dom- 
inates the expansion rate and the total perturbations. 
However, dark energy comes to dominate eventually, and 
then the potential starts to decay and the perturbations 
grow more slowly or start to decrease. This is also visible 
in Fig. [T] at very late times where the numerical solution 
for 6 starts to decrease. 

It is difficult to capture this behavior accurately. But 
in we introduced the variable Q(k,a) to describe the 
change of the gravitational potential due to the dark en- 
ergy perturbations. Q is defined through 



k 2 (f> = -A-n:Ga 2 Qp ri 



3aH, 



k 2 



-Vrr 



(19) 



If the dark energy or modification of gravity does not 
contribute to the gravitational potential (for example if 
the dark energy is a cosmological constant) then Q = 1. 
Otherwise Q will deviate from unity, and in general it is 
a function of both scale and time. 

Introducing the comoving density perturbation A = 
8 + 3aHV/k 2 , and looking at Eq. (T5|) we see that we can 
compute Q with 



Q-l = 



PdeA 



DE 



Pm A 7; 



(20) 



Just using the results during matter domination, we find 
that the resulting expression for Q is surprisingly accu- 
rate even at late times (see Fig. [5]). The reason is that 
both fluids, dark energy and matter, respond similarly 
to the change in the expansion rate so that most of the 
deviations cancel. We find that the sub-soundhorizon ex- 
pression below is accurate at the percent level, while on 
larger scales there are deviations of about 10 to 20% by 
today (depending on w). The latter can be corrected "by 
hand" in order to obtain a more precise formula, but the 
expressions are sufficiently accurate for our purposes and 
we keep them as they are. 

For matter, and during matter domination, the comov- 
ing density perturbation is extremely simple, A m = 5 a. 
For the dark energy we find 



A 



1 



DE 



1 + w 



" 1 — 3w 1 — 3w 



A., 



(21) 



on scales larger than the sound horizon. This means that 
the relative strength of the comoving density perturba- 
tions in the dark energy and the dark matter is constant 
on large scales, with those in the dark energy being sub- 
dominant for w < — for w close to — 1 the prefactor is 
approximately (1 +w)/4. Using the scaling of the energy 
density in matter and dark energy, we can derive that 



Q- 1 = 



1 - Sir, 



1 + <W 

1 - 3w 



, —3w 



ha~ 3w . (22) 



Here we defined a constant Qq since this expression will 
appear frequently later on. For f2 m = 0.25 and a — 1 Q 
interpolates smoothly between for w — — 1 and 3 for 
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w = 0. For w = -0.8 we have Q(a = 1) - 1 = Q « 0.18, 
that means that we do not expect more than about a 20% 
deviation of Q from 1 even on large scales, given current 
observational limits on w. We will keep using w = —0.8 
to illustrate what we can maximally expect to observe. 

On the other hand, on small scales the growth of 
dark energy perturbations is stopped by pressure sup- 
port. The dark energy perturbations stop growing once 
they are inside the sound horizon, while the matter per- 
turbations (for which there is no sound horizon) continue 
to grow. This leads to 



^DE 



2 V 



cx A m /a (23) 



on scales below the sound horizon. Here we neglected the 
dark energy velocity perturbations, since their contribu- 
tion to Ade is suppressed by a factor proportional to 
(Ha/k) 2 relative to Sue- From this expression, we then 
find that 



0-1 



l-o, 



3., .(Ha 
2 {1 + W) {-k 



-Sw 



(24) 



= (i-^la+w)^- 1 -*", 



where the result is more accurate when using the matter 
Ha as done above. Also on small scales, deviation of Q 
from 1 is maximal at late times for any dark energy that 
leads to acceleration, w < —1/3. But in general it will 
be suppressed relative to the large scales by the lack of 
growth of the dark energy perturbations once they are 
inside the sound horizon. 

For the variable Q we can construct a unified for- 
mula which accounts both for modes below and above 
the sound horizon, 



l-O, 



-3w 



(i 



l-3w 



2k 2 c 2 3 a 



(26) 



In Fig. [2] we compare Qtot with the numerical solution 
from CAMB for different values of the dark energy sound 
speed. The formula for Qtot interpolates between the two 
asymptotic regions. Close to the sound horizon it is not 
very accurate, but it is sufficient to use in e.g. Fisher- 
matrix codes to estimate the size of effects due to the 
presence of dark energy perturbations. 



V. IMPACT ON SOME OBSERVATIONAL 
QUANTITIES 

A. Growth and shape of the matter power 
spectrum 

What impact does the change in <p from Q ^ 1 have on 
the matter perturbations? Notice that the effect of the 
change in H from the onset of dark energy domination at 



<y io 




FIG. 2: The figure shows the behavior of the variable Qtot — 1 
for modes above and below the sound horizon. The solid lines 
are CAMB output (blue / upper line: <? a = 0.1, black / lower 
line: <? a = 1) and the dashed lines use Eq. (|26[) . red (upper) 
for Cg = 0.1 and yellow (lower) for <? s = 1. The vertical lines 
show the sound horizons, the left line for <? s = 1 and the right 
line for (? s = 0.1. The behavior on scales larger than the 
sound horizon is the same in both cases, but the growth of 
the perturbations slows down as the mode enters the sound 
horizon. As the modes enter the sound horizon earlier for 
(25) c s = 1 they stop growing earlier and stay smaller. We used 



again w — —0.8 and k = 200f/o 



late times is expected to be larger than the the effect from 
the presence or absence of dark energy perturbations, but 
here we want to quantify only the latter. We can write 
(j) — (j) m + 0de where <p m is the solution given by Eq. (fl"4"l) 
and 0de is (Q — l)0 m - Since the differential equation is 
linear, the solution with the total source is just the sum 
of the solution for each source (plus, as before, a decaying 
solution oc 1/a). The expressions we have found for Q 
are always power-laws in a, which guarantees a simple 
form for the velocity perturbation. 

Even though the dark matter does not have a sound 
horizon, the solutions will now depend on whether the 
k mode in question is larger or smaller than the sound 
horizon of the dark energy, simply because the deviation 
of <f> from (f> m depends on this. For fc-modes larger than 
the dark energy sound horizon, we find for the matter 
velocity perturbation 



V m 



S H 




(27) 



where w is the equation of state parameter of the dark 
energy fluid, and the constant Qq was defined in Eq. (|22p . 
At late times V m will deviate from this formula because of 
the change in the expansion rate, but we are again inter- 
ested in the impact of the perturbations, corresponding 
to the factor in curly brackets. The factor is less sensi- 
tive that the velocity itself, and it tracks the numerical 
result closely until dark energy domination sets in, giv- 
ing at least the right order of magnitude even today. For 
the optimistic case with w — —0.8 we find a deviation of 
about 6.5% at late times due to the presence of the dark 
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FIG. 3: The modification of the matter power spectrum P(k): 
The solid black line shows how the matter power spectrum is 
enhanced outside the sound horizon (vertical red dotted line). 
The blue dashed line shows the prediction from Eq. 1)280 for 
the modification outside the sound horizon. Here we used 



energy perturbations on scales where the matter pertur- 
bations become non-linear, especially if those scales are 
larger than the dark energy sound horizon. Second, for 
a dark energy model with a very large sound speed like 
Quintessence, for which c~ = 1, the matter power spec- 
trum will only be affected on the largest scales which have 
not been observed. It is obviously much easier to mea- 
sure the sound speed if it is low. However, the feature 
visible in Fig. [3] would provide a clear signature for the 
presence (and size) of a sound-horizon, especially if the 
corresponding change in the matter velocity field could 
be detected as well. 

A potentially more promising place to look for the im- 
pact of the dark energy perturbations is in way the mat- 
ter perturbations grow over time. Often the impact of 
the dark energy on the growth rate of the matter per- 
turbations is parametrised in terms of a growth index 7 
[HUlI, defined through 



-0.8 and c. 



10" 



dlog(A m ) 
dlog(a) 



fi(a) 



(29) 



energy perturbations. The fully numerical integration for 
this case gives about 4.5% change. The agreement is not 
as good as for the dark energy perturbations themselves, 
but still acceptable, especially since we are dealing with 
an indirect effect. 

This solution we can insert in the one for the density 
perturbation. In this case, </>' ^ 0, but since we are again 
dealing with a linear equation, we end up with a sum of 
three solutions, one for each source (the two parts of V m 
and <£>') as well as the constant peculiar solution: 



1 + ^Qoa~ 3w 



(28) 



-a 1 



Q a 



-3u> 



1 — 5w + 6w 



On large scales, larger than the causal horizon, the mat- 
ter perturbations are therefore enhanced by a factor 
1 + (3Qo)/2a -3 "\ which is of the same order as Q. On 
smaller scales, but still larger than the sound horizon, the 
factor is 1 + Qoa~ 3w / (I — 5w + 6w 2 ) which is smaller by 
a factor of 10 or so (depending on w). Inside the sound 
horizon, the growth of the dark energy perturbations is 
suppressed. 

It is customary to consider two different aspects of the 
matter perturbations: the power spectrum today and the 
growth rate of the perturbations. For the power spectrum 
today we can for example look at the ratio of 5^ with 
dark energy perturbations to 5^ without the perturba- 
tions. For our usual optimistic benchmark with w = —0.8 
we expect a 4% enhancement of P(k) on scales larger 
than the sound horizon. A numerical calculation shows 
that it is closer to 2%, see Fig. |3J 

Observationally there are two difficulties to detect the 
effect: first, it is unclear what happens with the dark 



We will now try to connect our expression (|28j) to 7 
through some rather crude approximations which should 
nonetheless be good enough to result in an order of mag- 
nitude estimate of the change of 7 due to the presence 
of the extra perturbations. On small scales (but always 
larger than the dark energy sound horizon) the differ- 
ences between S m and A m are suppressed by (Ha/k) 2 
and we have 



A m (fc » Ha) 



Using Al? = S a to model partially the late-time change 
of the expansion rate, and assuming that this term obeys 
the form of Eq. (f2"9"| with unperturbed growth index 70 
we find by performing the derivative 



5 a 



-3w 



1 — 5 1/7 + 6w 2 



(30) 



dlog(A m ) 
dlog(a) 



Aff'a 
(o) 



3wQqci 



-3tv 



A 



1 70 



Q a- 3w + (1 - 5w + 6uu 2 ) 
3(Q-1) 



5 — 6w — Q/w 



(31) 



where we used Eq. (f2"2"|) to reintroduce Q — 1. In general 
this cannot be cast in the form of m (a) 7 since the dark 
energy perturbations change the matter growth rate even 
during matter domination where f2 m (a) = 1. We can 
however connect to that form at least in the limit when 
^m( a ) just starts to deviate from unity. 

We would like to end up with SI (a) 71 with 71 = 70 + 
e. If e is small enough then we can use that f2 m (a) e « 
1 + elog(Sl TO (a)) and if we are close to Q(a) = 1 then 
additionally log(f2 m (a)) ~ Q m (a) — 1, giving finally 



e(fi m (a) - 1) 



3(Q-1) 
5 — 6w — Q/w 



(32) 
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According to [l8| (see also [Hj]), the growth index de- 
pends on Q through the combination 



7 

MQ) 



3(l-w~A{Q)) 
5 — 6w 

<9-i 



1 - Q m (a) 



so that in terms of our e above 

-3(Q-1) 



(33) 
(34) 

(35) 



which is sufficiently close to Eq. (|32[> . given the crudeness 
of the approximations used. Numerically the expression 
above is very close to Eq. (p?Tj) - closer than to the full 
numerical solution which remains more constant at late 
times. Both predict the correct deviation at early times, 
though. 

The denominator of the function A is 



(l-Q m (a))(5-6w) 



1 - Q m (a) = 



(1 - fi m )a 



-3w 



fl m + (1 — fl m )a' 



-3w 



(36) 



Not surprisingly, this is 1— O m today, so that A » (1+w), 
while at early times the second term in the denominator 
is suppressed so that for a <C 1 ^4 becomes Ade/A to m 
(1 + w)/4 on large scales. This has to be compared to 
1—w w 2. The dark energy perturbations thus change the 
growth index by a few percent or about 0.02 for w = —0.8 
on scales larger than the sound horizon. This will be chal- 
lenging to measure even by full-sky surveys like the pro- 
posed Euclid satellite mission which expects to achieve 
an error on 7 of less than this Q but only if enough 
modes can be measured, i.e. if the sound speed is close 
to zero. Once the perturbations enter the sound horizon, 
they stop to grow and so their impact on the matter per- 
turbations decreases and becomes rapidly negligible and 
thus impossible to detect. However, notice that this is 
purely the change of the growth rate due to the perturba- 
tions in the dark energy fluid! We defer a more detailed 
investigation of the detectability of the dark energy per- 
turbations in the dark matter power spectrum to a later 
publication 



B. The integrated Sachs- Wolfe effect 

After the matter power spectrum, we focus on the inte- 
grated Sachs- Wolfe (ISW) effect for cosmological models 
with non-zero contribution from the dark energy pertur- 
bations. The motivation is that the ISW part of the 
spectrum is the most affected by the dark energy, see 

The ISW effect results from the late-time decay of 
gravitational potentials. The total blueshifting or red- 
shifting of the CMB photons caused by the change in 
the potential during the passage of the photons induces 
temperature fluctuations [24| : 



T J Or j 



a'H^dx, (37) 
oa 



where r denotes conformal time. In the last step, we 
have replaced the integration variable by the comoving 
distance x which is related to the conformal time by dx = 
—cdr = —cdt/a; here we assume a zero anisotropic stress 
component for all the species in the Universe. 

In Fourier space, the derivative of the gravitational po- 
tential with respect to the scale factor a can be expressed 
as: 



3 H Cl r , 



2 ak 2 



{Q(a,k) A' m (a,k) + 



+ Q' (a, k) A m (a, k) - -Q (a, k) A m (a, fc) j . (38) 

where the prime denotes the derivative with respect the 
scale factor. There are three main ways in which the 
dark energy perturbations change the ISW effect: They 
change A m itself, as discussed in the last section, both 
the shape (a very small effect) and the growth. They 
change cp' additionally through the presence of Q in the 
last and first term of Eq. (|38[) and their variation enters 
as well through Q'. 

In linear perturbation theory all k modes evolve inde- 
pendently, so that we can decompose the dark matter 
density contrast as 



A m (a, k) = aG(a, k)A m (k). 



(39) 



Here A m (fc) = A m (a — l,k). The so-called growth 
factor G is usually independent of scale: during mat- 
ter domination G is constant and the late-time change 
in the expansion rate affects all scales equally. However, 
the contribution from the dark energy perturbations in- 
duces a scale dependence because of the existence of the 
dark energy sound horizon, see for instance (25j . For the 
growth rate, we will use the 7(0, fc) from Eqs. (|33[) and 
(I34p . which provides a sufficiently good approximation. 
We can then write Eq. (|3"5)) as: 



3 i?o^m d 



2 fc 2 da 



{G(a,k)Q(a,k)^A m (k). (40) 



The line of sight integral for the ISW-temperature per- 
turbation £ can now be written as 



c 



d X W c ( X ) A m (fc) 



(41) 



where we introduced the weight function: 



3 H n Q r , 



H^-{G(a,k)Q{a,k)} (42) 



which allows the expressions for the ISW-auto spectrum 
CcC W t° be written in a compact notation, applying the 
Limber-projection [26j in the flat-sky approximation, for 
simplicity: 



dx- 



W? (x) 



c 

X 2 



Pa a (fc = t/x)- 



(43) 
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FIG. 4: ISW power spectrum for c? a — 10 4 (red solid line) 
and Cg = 1 (blue dashed line) for a model with w = —0.8. 



Paa (k) is the linear matter power spectrum today, 
which can be written as: 



fc 3 P A A (ji) 

2tt 2 



S 2 h 



(44) 



Here 8h is the amplitude of the present-day density fluc- 
tuations at the Hubble scale and T (k) is the transfer- 
function for CDM. We adopt the fit by Eisenstein & Hu 
J27[. We have neglected the additional impact from the 
dark energy perturbations onto the dark matter power 
spectrum as it is only of the order of a few percent (see 
the discussion in the last section). The Eisenstein & Hu 
fit agrees with the CAMB output to a precision of about 
4% for both low and high values of c 2 , more than suffi- 
cient for the purposes of this section. 

In Fig. H we plot £(£+ 1) C £ /{2n) using Eq. (g3]) for 
two different values of the dark energy sound speed. The 
differences between the two curves come from the term 
(d(QG)/da) 2 . Since the ISW power spectrum depends 
on the derivatives of the product of the growth factor 
G(a, k) and Q(a, k) we need to look at 



(GQ)' = G'Q + GQ' 



(45) 



The deviation of Q from 1 is never enough to explain the 
differences between the two curves in Fig.|4j Indeed, tak- 
ing Q — 1 while keeping Q' barely changes the results. 
The relative size between the first and second term de- 
pends mainly on Q'G/G' . From Eq. (f2"6"|) we see that 
aQ' oc (Q — 1) with a proportionality factor of about 
2 to 3. This is a small number, but it is boosted by 
G' < G, since aG'/G = fl m (a)^ - 1. The dark energy 
slows down the growth of the dark matter perturbations 
so that G' < 0. However, Q' > because of the rel- 
ative increase of the dark energy density enhancing the 
importance of the dark energy fluctuations. The two con- 
tributions will partially cancel and so decrease the result, 
unless one term dominates strongly, e.g. for modes inside 
the sound horizon. This is the reason why the ISW con- 
tribution to the CMB power spectrum decreases as the 
sound speed increases. 



FIG. 5: 

and w 



We plot the magnification factor A 2 with k = 200-ffo 
= —0.8 for two different values of the sound speed: 



10 (red solid line) and c s = 1 (blue dashed line). The 



vertical line at a — 0.62 shows when the A;-mode enters the 
sound horizon for (? s — 10 ~ 4 
mode enters at a — 0.62 x 10" 



(for the high sound speed the 
- 4 )- 



We can illustrate the effect by looking at the gravito- 
power spectrum, 



P<b>&> — 



\ 2k 2 



) 2 - 

G{a,k)Q(a,k)\ P AA (fe) . 



da 



(46) 

If we assume that the dark energy does not cluster (i.e. 
Q = 1), Eq. (H| reads: 



Prh'&t — 



( 3Hq O m \ f d 



V 2k 2 J Paa(*). (47) 



and the growth factor becomes again a function of time 
only. 

Comparing Eq. P6"]) and Eq. HT)) we can define a mag- 
nification parameter for the ISW power spectrum: 



A 2 = 



d(G(a,k)Q(a,k)) /da 
dG(a)/da 



(48) 



This is precisely the expression discussed above, A = 
Q + Q'G/G 1 . 

In Fig. [5] we plot the magnification factor for two val- 
ues of c 2 for k = 200ffo. For a sound speed equal to 1 the 
dark energy perturbations enter the sound horizon very 
early and stay small until today, even when taking into 
account that the dark energy density grows relative to 
the dark matter density. As expected, they do not affect 
the ISW effect significantly in this case. For small sound 
speeds the dark energy perturbations partially cancel the 
contribution from G' and decrease A by about 30% for 
scales above the dark energy sound horizon. Their slower 
growth inside the sound horizon leads to a smaller can- 
cellation, as can be seen in Fig. \5\ Based on these obser- 
vations, we expect that a low sound speed decreases the 
ISW effect (which goes like A 2 ) by about 50%, consistent 
with Fig.|U Overall, it seems that Q 1 can provide a more 
sensitive probe of the dark energy perturbations than Q. 
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VI. CONCLUSIONS 

In this paper we have derived simple analytical expres- 
sions for the behavior of the perturbations in fluid dark 
energy models with vanishing anisotropic stress, equa- 
tions (fT5|) to (fT8"|) . These expressions are valid for models 
with constant w < and cj, > 0. They were derived 
under the assumption of matter domination, but they 
allow to compute the function Q{k,a]w,c 2 s ) which de- 
scribes the deviation of the Poisson equation from the 
case without dark energy perturbations and which is rel- 
atively insensitive to the late-time deviations from matter 
domination. The expressions for Q are given in Eq. (|22p 
for scales larger than the sound horizon and Eq. (|2"5|) for 
scales smaller than the sound horizon. We also give a 
single interpolating equation in Eq. (|26p which is useful 
for Fisher-matrix calculations that include dark energy 
perturbations. In models without anisotropic stress Q 
completely characterizes the dark energy perturbations 
and represents a fingerprint that allows to differentiate 
between different models with the same background ex- 
pansion rate but a different evolution of the perturba- 
tions. 

We expect our results to hold generically for scalar- 
field like models. Large changes can appear because 
of rapidly varying (especially oscillating) w leading to 
resonance-like behavior (see e.g. the rather contrived 
"phaxion" model in [Hj]), for models with non-zero effec- 
tive anisotropic stresses (like DGP 0,[11|) or couplings 
between the dark energy and the dark matter [2II l3fl l31j| . 
so such models will lead to different results. 

We find that the dark energy perturbations are always 
smaller than the perturbations in the dark matter, at 
least by a factor (1 + 10), but on scales larger than the 
dark energy sound horizon they are only suppressed by 
an additional factor of order unity. The dark energy per- 
turbations do not grow on scales smaller than the sound 
horizon, so that these perturbations are much more sup- 
pressed relative to the dark matter. The impact of the 
dark energy perturbations on the gravitational potential 
is additionally influenced by the relative energy density, 
so that measurable deviations tend to appear only at late 
times. 

To demonstrate the usefulness of the equations we then 
used the formula for Q to investigate the change in the 
dark matter power spectrum and the ISW effect if w = 
—0.8. The changes in the matter power spectrum at late 
times are of the order (Q(a, k) — l)/5 which corresponds 
to a few percent on scales larger than the sound horizon 
of the dark energy. The growth index 7 is changed by 
about 0.02 on the same scales. We further find that the 
impact on the ISW effect is due to the growth of the 
dark energy perturbations, Q' having a larger effect on 
<f>' than naively expected from the size of Q. It is much 
larger than the impact of the dark energy perturbations 
on the matter power spectrum, but because of cosmic 
variance it is more difficult to constrain observationally. 

The dark energy perturbations of the class of mod- 



els investigated here vanish as w — > — 1. However, if w 
is different from —1 and especially if the sound speed 
of the dark energy is less than the speed of light, then 
there is hope that the effects from the dark energy per- 
turbations could be seen with cosmological observations. 
Although very challenging, it is nonetheless worth the 
effort since the perturbations provide a much more pre- 
cise "fingerprint" of the dark energy than the equation 
of state parameter w. 
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APPENDIX A: COMPARING WITH CAMB 

Since we are dealing with linear perturbation theory, 
all fc-modes evolve independently and for each fc-mode 
the perturbations depend linearly on the normalization 
given by the constant <5 . As we can choose this constant 
arbitrarily for each fc-mode, we should really think of it 
as S (k) . In the standard cosmology, its value is set in the 
very early universe by inflation. This scenario then pro- 
vides the initial conditions for the differential equations. 
However, when we compare our results with numerical 
solutions from e.g. CAMB, we have to take into account 
as well that the Universe was radiation dominated at 
early times. Since both the expansion rate of the Uni- 
verse and the dominant contribution to the gravitational 
potentials are different during radiation domination, we 
expect to find a different behavior for the matter and 
dark energy perturbations. It turns out that the change 
in the expansion rate strongly suppresses the growth rate 
of the perturbations, leading to an only logarithmically 
growing solution (32j. 

This is not directly relevant to our solutions as we are 
content to limit our expressions to the matter dominated 
and later epochs. One nuisance is that we have to disre- 
gard radiation domination in some of the figures where we 
prefer a larger dynamical range to show the (formal) evo- 
lution of the perturbations in the different regimes. We 
note this in the figure captions where applicable. Another 
issue concerns the normalization for comparison with nu- 
merical codes: as the perturbation growth of sub-horizon 
modes is delayed by the radiation dominated epoch, the 
k modes which enter the horizon during that period end 
up with a lower normalization than expected if they are 
normalised at early times. The pragmatic solution here 
is to normalize the perturbations instead in the late uni- 
verse, after the onset of matter domination, 

Sm (a = a 1 ) = d in = 5 L+^^j (Al) 

which fixes S (k) in terms of 5i n {k) at a given scale factor 
ai . We also note that radiation pushes the growth of both 
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dark matter and dark energy perturbations to later times 
so that the ratio (which is relevant for Q) is basically 
unchanged. 

If a more detailed treatment is desirable, then this can 
be obtained by following the discussion in standard texts 
like e.g. (32J- We prefer to avoid these additional com- 
plexities here since they are not required for the main 
points of our work. 

In this paper we are working in the Newtonian gauge, 
but CAMB uses the synchronous gauge. In order 
to compare the CAMB output with our formulae, we 
need to transform the relevant quantities. The energy- 
momentum tensor Tjf (Syn) in the synchronous gauge 
is related to the (Con) in the conformal Newtonian 
gauge by the transformation: 

T " {Syn) = d^d¥ T ^ C(m) (A2) 

where x 11 and x^ denote the synchronous and the con- 
formal Newtonian coordinates respectively. The relevant 
transformations then are II ill 



6 (Syn) 

V(Syn) 
5p (Syn) 
a (Syn) 



8 (Con) — ar: 
P 

V (Con) - a(l + w)k 2 
Sp (Con) — oip 
a (Con) 



(A3) 

(A4) 
(A5) 
(A6) 



where p is the energy density at background, 8 is the 
density contrast, V is the velocity perturbation, Sp is 
the pressure perturbation, a is the anisotropic stress and 

the function a = x° — x° = (h + i)j /2k 2 . Here n and 

h are two scalar fields characterizing the scalar modes 
of the metric perturbation in synchronous gauge. This 
transformation also applies to individual species when 
more than one particle species contributes to the energy- 
momentum tensor. 



APPENDIX B: DECAYING MODES 

In the main text we were able to find simple approxi- 
mate solutions to the full system of equations by neglect- 
ing certain terms, especially those linking S and V for 
the dark energy perturbations. This works very well for 
the parameter values of interest to us, specifically w < 
and < c 2 < 1. For c 2 s < we expect rapid growth of 
the density perturbations from a well-known instability. 
This instability, however, requires exactly the coupling 
which we neglected. In this appendix we want to take a 
closer look at the "decaying modes" and their behavior 
especially for c 2 < 0. For this purpose we need to revert 
to the full system. To allow for consistent simplifications 
we can then cast it in the form of a single second-order 
differential equation by combining Eqs. ^ and ([7]): 



6" 



-(1-w) 
a 



IT 



A[_ 

T 



5' 



(Bl) 



Ak 2 c 2 
(a 2 H) 2 



(2-3c 2 s ) + — -a- 



A' 



A 



(l + w)A 2 

9— k 

(a 2 H) 2 







where: 



.4 



A' = 



9a 2 H 2 



k 2 

9 (2aH 2 + 2a 2 HH') (c 2 
k 2 



(B2) 
(B3) 



We expect to see rapid perturbation growth for imagi- 
nary sound speeds on small scales. For k 3> 1 then A ~ 1 
and A' ~ 0; Eq. (|BT|) becomes: 



5" + — (l-2w) 6' 
2a 



+ 



— (c 2 

2a 2 { s 



(a 2 H) 2 



(B4) 



0. 



At the sound horizon c 2 k 2 = (aH) 2 = H^Vlm/a so that 
for sub-sound-horizon modes v(a) 2 = k 2 c 2 a / (H^Vl m ) ^> 
\(c 2 — w)(l — 6c^)|, and we can neglect the contribution 
from the second term, v quantifies how deep the mode 
is inside the sound horizon, and it will become complex 
for c 2 < (notice also that v grows as ^fa). We are then 
left with: 



I' 1 



2w) 6' 



k 2 c 2 



(a 2 H) 2 



(B5) 



The full solution to this equation contains the one found 
earlier, Eq. (fT?)) . as well as two additional ones, 

<5i cx v(a)- n J n (2v(a))T(Z/2-Zw) (B6) 
5 2 oc is(a)- n J- n (2v(a))T(l/2 + 3w) (B7) 

where n = (1 — 6w)/2. If c 2 < then the absolute value 
of the Bessel functions will grow exponentially fast since 
its argument is complex, as expected for the instability. 
For super-horizon modes (fc <C 1) Eq. (|B2[) reduces to 
— w) jk 2 and Eq. (|B1|) becomes: 



8"- 



5 - &w gl , 9 (c 2 - w) x | 9 (1 + w) (c 2 - w) 



2a 



2a 2 



i 2 k 2 



k 2 4> = 0. 
(B8) 

The full solution to this equation contains the one 
found earlier, Eq. (fTB)) . as well as the following ones, 



(B9) 
(B10) 



x § (-l+2w-y/-8c 2 s + (l+2w) 2 ) 

01 (X a 4 V V 'J 

„ f (-1+2W+J -Sc 2 3 + (l+2w)A 

02 oc a 4 V v J 
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The interesting term here is the square root, in order 
to find a growing solution it needs to be real: —8c 2 . + (1 + 
2w ) 2 > 0. This is the case for 

(l + 2w) 2 >8c 2 . (Bll) 



This means that if we want to keep the sound speed pos- 
itive we need to have an equation of state parameter that 
is either very negative or well larger than —1/2. In either 
case the solutions are decaying for w < 0. 
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